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Two-Component Nonlinear Schrodinger Systems 

Tai-Chia Lin *^and Juncheng Wei * 



Abstract 

Recently, skyrmions with integer topological charges have been observed numer- 
ically but have not yet been shown rigorously on two-component systems of nonlinear 
Schrodinger equations (NLSE) describing a binary mixture of Bose-Einstein condensates 
(cf. [2] and [25] ). Besides, half- skyrmions characterized by half-integer topological charges 
can also be found in the nonlinear a model which is a model of the Bose-Einstein con- 
densate of the Schwinger bosons (cf. [IS])- Here we prove rigorously the existence of 
half- skyrmions which may come from a new type of soliton solutions called spike-vortex 
solutions of two-component systems of NLSE on the entire plane M. 2 . These spike- vortex 
solutions having spikes in one component and a vortex in the other component may form 
half-skyrmions. By Liapunov- Schmidt reduction process, we may find spike-vortex solu- 
tions of two-component systems of NLSE. 

1 Introduction 

Spikes and vortices are important phenomena in one-component nonlinear Schrodinger 
equations (NLSE) having applications in many physical problems, especially in Bose-Einstein 
condensation. In the last two decades, there have been many analytical works on both spikes 
and vortices, respectively. One may refer to pj5] for a good survey on spikes, and p], [H] and 
[21] for survey on vortices. Recently, a double condensate i.e. a binary mixture of Bose-Einstein 
condensates in two different hyperfine states has been observed and described by two-component 
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systems of NLSE (cf. [22]). It would be possible to find spike and vortex solutions from two- 

component systems of NLSE. However, until now, there is no result to deal with spike-vortex 
solutions having spikes in one component and a vortex in the other component. In this paper, 
we want to find such solutions and investigate the interaction of spikes and vortices. 

Ordinary skyrmions being of topological solitons resemble polyhedral shells which look like 
closed loops, possibly linked or knotted. Besides their intrinsic fundamental interest, skyrmions 
have important applications in nuclear physics (cf. p2]), and analogous structures are postu- 
lated for early universe cosmology (cf. [1]). To get a skyrmion, one may consider the multi- 
component wave function which may introduce the extra internal degrees of freedom and result 
in a nontrivial structure characterized by topological charges. For two-dimensional Skyrme 
model, skyrmions have been investigated by the method of concentration-compactness (cf. [IB]). 
In a double condensate, skyrmions with integer topological charges have been observed by 
numerical simulations on two-component systems of NLSE (cf. [2] and [25]). Besides, half- 
skyrmions characterized by half-integer topological charges can also be found in the nonlinear 
a model which may describe the Bose-Einstein condensate of the Schwinger bosons (cf. [T8]). 
Here we prove rigorously the existence of half- skyrmions in a double condensate using spike- 
vortex solutions of two-component systems of NLSE. 

To get spike-vortex solutions, we study soliton solutions of time-dependent two-component 
systems of NLSE as follows: 



where the spatial dimension n = 2 and ipj = ipj(x,t) G C for j = 1,2. The system (11. ip 
is a standard model to describe a double condensate. Physically, ip/s are the corresponding 
complex-valued wave functions, and the coefficients ~ — for i,j = 1,2, where a^'s and 
O12 = 0,21 ar e the intraspecies and interspecies scattering lengths, respectively. When the spatial 
dimension is one, i.e. n = 1, it is well-known that the system (II. ip is integrable, and there 
are many analytical and numerical results on soliton solutions of coupled nonlinear Schrodinger 
equations (e.g. [7], [8], [9]). Recently, from physical experiments (cf. [3]), even three-dimensional 
solitons have been observed in Bose-Einstein condensates. It is natural to believe that there 
are two-dimensional (i.e. n = 2) solitons in double condensates. However, when the spatial 
dimension is two, the system (II. ip becomes non-integrable and has only few results on two- 
dimensional solitons. This may lead us to study two-dimensional soliton solutions of the system 
(II. ip and find different types of solitons. 

In the vicinity of a Feshbach resonance, scattering lengths a^'s depend sensitively on the 
magnitude of an externally applied magnetic field (cf. [23] and [21]), allowing the magnitude and 
sign of flifs to be tuned to any value. Generically, when both Pj/s are positive, the system (11. ip 
is of self-focusing and has bright soliton solutions on the associated two components. On the 
other hand, when both /3j/s are negative, the system (II. ip is of self-defocusing and has dark 
soliton solutions on the associated two components. Here we have interest on the case that (3n 
and @22 have opposite signs which may result in a new type of soliton solutions called spike- 
vortex solutions of the system (II. ip i.e. one component has spikes and the other component 
has a vortex. Furthermore, we may obtain half- skyrmions by these spike-vortex solutions. 
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To obtain soliton solutions of the system (11.11) . we set 

ipj{x,t) = e^ 1 ^ ■ Uj{x) , uj 6 C, j = 1,2. (1.2) 

Substituting (11. 2ft into ( 11.11) . we may obtain a two-component system of semilinear elliptic 
partial differential equations given by 

Aui - Aiui + f3 n |ui| 2 ui + (3i2 ui \u 2 \ 2 = 0, ^ 
Au 2 - X 2 u 2 + P 22 \u 2 \ 2 u 2 + (3i 2 \ui\ 2 u 2 = 0, 

where fii 2 is the coupling constant. In [IS]- [IS], we studied the ground state solutions of (II. 3p 
for the case that fin > 0, [3 22 > and u/s are positive functions. Namely, we investigated the 
following problem: 

Au — Xiu + /SiiM 3 + Pi 2 uv 2 = 0, 
Av - X 2 v + (3 22 v 3 + (3 12 u 2 v = 0, 
u,v > 0, u,v E H 1 '^ 

Due to each Xj > and /3jj > 0, both u and v components have attractive self- interaction 
which may let spikes occur in these two components. One the other hand, when n = 2, Xj < 
and (3jj < 0, j = 1, 2, and u/s are complex-valued solutions of (II. 3p . vortices may exist in both 
Mi and u 2 components (cf. [T2]). 

In this paper, we study the case that n = 2, Ai,/5n > 0, A 2 , /3 22 < 0, and Ui is positive but 
«2 is complex-valued function. Without loss of generality, we may assume that 

Ai = — A2 = /?n = — j3 22 = 1 . 

Namely, we study the following system: 

Au - u + u 3 + (3u\v | 2 = inlR 2 , 
Av + v- \v\ 2 v + pu 2 v = in M 2 , 

where n > and v G C. To get skyrmions, a defining property of the skyrmion is that 
the atomic field approaches a constant value at spatial infinity (cf. [26]). Hence we may set 
boundary conditions of the system (11.51) as follows: u(x) — > and \v(x)\ — > 1 as |x| — » 00. 

As = 0, the first equation of the system (11.51) becomes a standard nonlinear Schrodinger 
equation given by 

A«-m + m 3 = 0, ueH\R 2 ) (1.6) 

which has a unique least-energy spike solution w = w(r),r = \x\ satisfying w'(r) < for r > 
and 



w(r) = A n r 2e r 



w'(r) = — A r 2 6 r 



;i-7) 



Actually, this is also a typical spike solution of nonlinear Schrodinger equations. On the other 
hand, as (3 = 0, the second equation of the system (11.51) can be written as 

Av + v- \v\ 2 v = 0, v = v(z) G C for zeC = l 2 , (1.8) 
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which is of conventional Ginzburg-Landau equations (cf. PQ) having a symmetric vortex solution 
of degree d G Z\{0} with the following form 

v d (z) = S d (r)e^ M , (1.9) 

where S d (r) satisfies 

I + \S' d - fjSd + S d - S d = 0, 



S d (0) = 0, S d (+oo) 



and 



S' d (r)>0, 5 d (r) = l-^ + 0(i) asr-,+oo. (1.11) 

Here we want to prove that when /3 increases or decreases, there exist spike-vortex solutions of 
the system (ll.5p . This may become the first paper to illustrate such solutions of two-component 
nonlinear Schrodinger systems. 

The main purpose of this paper is to construct a spike- vortex solution (u,v) of the sys- 
tem (ll.5p such that u ~ w and v ~ v d . Actually, the main difficulty of this paper is to study 
the interaction of typical spike and vortex solutions. Our first result shows that even a weak 
repulsive interaction (j3 > being small) can produce abundant bound states by solving the 
system (ll.5p . More precisely, we have 

Theorem 1.1. Let n = 2, d g N and k > 2 s«fe/y 

(%) k > 2 is any positive integer if d — 1, 

(nj 2(d — 1) ^ mod k i.e. there does not exist any integer fi such that 2(d — 1) = kfi if 
d>2. 

Then for (3 > sufficiently small, the problem M.5\) has a solution {up, Vp) satisfying up(z) > 
for z G C, vp(0) = with degree d, and up(z) — > 0, —> 1 as \z\ — > oo. Moreover, as 

(3 — > 0+, (up,vp) has the following asymptotic form 

k 

up(z) = J2w(z-&+0(\(3\), (ii2) 
vp(z) = S d (r)e^ lde e^^ z \ ^p{z) = 0(\(3\) G C, 
where w is the unique radial solution of $1.6\) . (£f, ■ forms a regular k-polygon with 



n i — rSjrfjf— 1) 

^ = i/se V=r-^ | j = i, (1.13) 
and — > +oo as P ^ 0+ satisfying lp = lp + 0(1), where lp satisfies 

A picture of (tig, i^) is given by 
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In [15]- [16], the positive sign of (3 may give inter-component attraction when the inter- 
component interaction is only for spikes. Conversely, the positive sign of (3 may contribute 
inter-component repulsion when the inter-component interaction is for spikes and a vortex. 
The inter-component repulsion between u and v components may balance with self-attraction in 
u-component so a new kind of soliton solutions called spike-vortex solutions of two-component 
nonlinear Schrodinger systems can be obtained in Theorem 11.11 

For (3 < 0, we may consider the radial solution of (11.51) given by 

u = u(r ) , v = f{r) e^™ , (1.15) 

satisfying 

u" + i u' - u + u 3 + (3f 2 u = , Vr > , 

f" + lf'-^f + f-f 3 + (3u 2 f =0, Vr>0, . . 

u'(0) = 0,m(+oo) = 0, 1 ' 

/(0) = 0,/(+oo) = l, 

where d E N and (r,6) is the polar coordinates in M 2 . Then we have the following existence 
theorem: 

Theorem 1.2. A ssume that n = 2 and (3 < 0. Then the problem U.5\) has radially 
symmetric solutions of the following form: 

u = u{r), v = f{r)e V ^ Tde (1.17) 

where u(r) is strictly decreasing, f(r) is strictly increasing and (u, f) satisfies $1.16}) . 

Remarks: 

1. It is easy to see that the solution (w(r), /(r)) is unique for the system (I1.16P if (3 is small 
enough. It is an interesting question to study the uniqueness for general (3 < 0. 

2. Note that the solution (w(r), f{r)e"^ de ) is not a global minimizer for the corresponding 
energy functional of (I1.5P since the equation of u is superlinear. It is conjectured that there 
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exists a least energy solution (see definitions in [IE] ) satisfying E[u(r), f{r)e^^ de ] < E[u,v] 
for any solution (w, v) of (jl.5p with deg(w) = d. Here the corresponding energy functional is 
defined as 

1 /" /m |2 i „,2\ 1 / „,4 , 1 / |t-7„.|2 , ^ /" /i l„.|2\2 ^ /" „,2|„.|2 



E[u,v] = - (|V«r + ^)- z / u 4 + -/ |V^ + - (l-MT-77 (1-18) 
To get half-skyrmions, we may define a S" 2 -valued map 




\/ M 2 + v\ + 

where (u, 1 ^2) is the spike- vortex solution obtained in Theorem 1 1 . 1 1 and \T72[ Generically, 

the topological charge of S 2 - valued maps is defined by (cf. [6]) 

Q — — I ~n ■ (d x ~n A d y Ti) dx dy . (1.20) 



47T _ 

Then we have 

Theorem 1.3. TheS 2 -map defined by U.19\) is of half-skyrmions with topological charge |. 
Throughout the rest of the paper, we assume that 

h--y<l<h + l (1.21) 

for some suitable 7. Note that 

where c& is constant depending on k only. Besides, unless otherwise stated, the letter C will 
always denote various generic constants which are independent of j3, especially for j3 sufficiently 
small. The constant a G (0, |) is a fixed small constant. 

The rest of this paper is organized as follows: In Section 2, we introduce useful properties 
about the spike solution w and the symmetric vortex solution Vd- In Section 3, we define the 
approximate solutions of spike-vortex solutions and derive the associated error estimates. In 
Section 4, we use Liapunov-Schmidt reduction process to find spike-vortex solutions. Then we 
may complete the proof of Theorem 11.11 and 11.21 in Section 5 and 6, respectively. Finally, we 
give the proof of Theorem 11.31 in Section 7. 

Acknowledgments: The research of the first author is partially supported by a research Grant 
from NSC of Taiwan. The research of the second author is partially supported by an Earmarked 
Grant from RGC of Hong Kong. 
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2 Properties of Spikes and Vortex 

We recall some properties of w and Sd(r)e y ^ de , where |d| = 1 or \d\ > 1 and 2(d —1)^0 
mod k. Let 



Li[0] = A0-0 + 3w 2 0, 

L 2 [ip] = Aip + -0 - Sfr() - 2Re(S d (r)e~^ Tde ip)S d (r)e^™ 



(2.1) 



for is a real-valued function and ip is a complex-valued function. For convenience, we may 
define the conjugate operator of L 2 by 



-ue 



L 2 := e-v-^L^ 
Then by simple computations, it is easy to check that 

^2 [^l + v 7 ^'^] = £2,1 [^1, ^2] + ^/-L L 2i2 [ipi, 2 ] , 
for ^1 and 0> 2 are real-valued functions, where 

d 2 2d 

L 2jl [V>i, 2 ] = A^i + (1 - 35j)^i - -^1 - -^2 



(2.2) 



(2.3) 



d 2 2d 

£2,2^1, $2] = A^ 2 + (1 - S 2 )</; 2 - — ^ 2 + —d e il) 



1 ■ 



Set a function space 



-1 ■ 



<t>(z), 



z 



(p{ze y k 

^{ze^ 1 ^) = e^^ip^z), tp(z) = tp(z)* 



(2.4) 



where k > 2 is an integer. Hereafter, both the over-bar and asterisk denote complex conjugate. 
We remark that the equation (11.51) is invariant under the following two maps 



x ze^ 











_ -1 Z7T 
1 L. 



* ^(ze^ 



(2.5) 



Therefore, we may look for solutions of (11.51) in the space J^. We first have 

Lemma 2.1. 

(1) Suppose Lx[0] = 0, G # 2 (M 2 ) and 0(z) = 0(z). T/ien 0(z) = c§^{z), where z 



z\ + v — 1-^2, 2j G R c is a constant. 
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(2) Suppose 

L 2 [tfj] = 0, |V| < C, ip(z)* = tp(z), and ^(ze^ 1 ^) = e^^^(z) , (2.6) 

where C is a positive constant. Then ip = 0, provided that \d\ — 1 or \d\ > 1 and 
2(d- 1)^0 mod fc. 

PROOF. (1) is easy to show. See Appendix C of [20]. We only need to show (2). We 
firstly state the proof for the case when \d\ = 1. For simplicity, we may assume d = 1. 
From |21j (Theorem 3.2), we know that 

^ = c (V^TS(r)e^ 6 ) + A( 5(r ) e V=») , (2.7) 



dz 3 

^0,0 



^0,7 



where Cj's are constants. It is easy to calculate that {^/ —lS^e^^^j = —y/—lS(r)e v ^ 6 ' ) 

lfoi= (^)'£ + ^ + V=T(^)'*V, and ^o, 2 = (^)'f + v^T[(^)V + ^], where x = z x 
and y = z 2 . Consequently, (-0 o ,o(^))* = -il>o,o{z), (ipo tl (z))* = ^ ,i(z) and (^0,2(2))* ^ 4>o,2( z )- 
Moreover, due to ip(z)* = ip(z), we have cq = c 2 = 0. Hence we only have 

if)(z) = cii> 0>1 (z) . 

However, it is obvious that ^ ,i doesn't satisfy ^ ,i (ze^^^) = e^~^^ ' ^0,1(2) ■ Thus C\ = 
and ijj = 0. 

Now we give the proof for the case that |d| > 1. From [TJ], the solution ip locally may 
become a linear combination of ipd,o(z) = h(r) e^~^ de and the following forms 

ipd,m( z ) = a ( r ) e^ 1 ^-^ 6 + b(r) e V=T(d+m)e ; (2.8) 

for m G N, where z = r e^~^ B and (r, #) is the polar coordinate. Here h, a and b are real-valued 
functions. Actually, these forms are invariant to the operator L 2 so one may decompose the func- 
tion space £ as invariant subspaces having the forms like if)a,o and ipd,m^- Then the condition 
ip(z)* = ip(z) may imply ipd,o(z)* = ipd,o(z)- However, since ipd,o( z ) — yf—lh(r)e yCZ ^ de , then 
ipdfiiz)* — —ipd,o(z) which gives ipd,o(z) = 0. Besides, the condition ip^ze"^^) = e^^^~i/j(z) 
may give ^d^ze^ 1 ^) = e^T^^z). Consequently, 



Hence a = or b = if e V=T(d-m-i)^ ^ 1 qt e ^i [d+m -i)^ ^ L Due tQ £ 2 [^ m ] = , a (r) 
and 6(r) satisfy 

a" + \a'-^a + {l-Sl)a-{a + b)Sl =0, 
6" + i6 , -^6+(l-53)6-(a + 6)S3 =0. K ' ) 



This implies a = & = 0ifa = 0or& = 0. Thus a = b = if e V=i(d-m-i)% ^ j or 

e V=T(<i+m-l)^ _^ L J t ig triyial that e J=Hd-m-l)% ^ 1 Qr e V=T(d+m-l)^ ^ 1 if 2(d - 1) ^ 

mod /c. Therefore a = 6 = i.e. ipd,m = if 2(d — 1) ^ mod k, and we may complete the 
proof of Lemma 12.11 □ 

To study the properties of L 2 (or L 2 ), we introduce some Sobolev spaces. Let a G (0, |). 
We introduce Hilbert spaces X a and Y a as follows: 



(i + |x| 2+Q )(M + H) 2 <+oo 



u = u 1 + V^Tu 2 e Ll c (R 2 ;C) 
equipped with inner product (u,v) Xa = (1 + \x\ 2+a )(uiVi + u 2 v 2 )dx, and 



Vl 



-lv 2 e W, 



2,2, 
loc \ 



2 ;C) 



|Av| 2 (l + |x| 2+Q )da; + 





V 


2 


1 + 


\x 


2+a 



equipped with inner product (u, v)y a 
to the inequality 

\h\ < 



(Au,Av) Xa + 



u ■ V 



\x 



2+a 



dx < +oo 



dx, respectively. Thanks 



(i + M) 2+c 

Besides, we see that X a has a compact embedding to L 1 (R 2 ). Originally, these spaces are real- 
valued function spaces introduced in Chae and Imanuvilor (cf. [S]). Here we generalize X a , Y a 
as complex-valued function spaces, and regard the operator L 2 from the space Y a to the space 
X a . We list some properties of X a and Y a , whose proofs are exactly the same as in [5]. 

Lemma 2.2. 

(1) Let v G Y a be a harmonic function. Then v =constant. 

(2) \/v G Y a , we have \v(z)\ < C x j |^| (ln(l + \z\) + 1), \/z G 1R 2 . 

(3) The image of L 2 (or L 2 ) is closed in X a n S , where S = {ip = 4>{z) G C : (0, ip) G S} . 

Now we study the invertibility of L 2 (or L 2 ) on the space Y a PI S . 

LEMMA 2.3. For a G (0, |). Then operator L 2 (or L 2 ) from the space Y a n So onto the 
space X a fl Sq is invertible. Furthermore, 



Ya 



<c\\L 2 mi 



<C\\L2M: 



(2.1i; 



PROOF. It suffices to consider the invertibility of L 2 . Then the invertibility of L 2 follows 
from a trivial transformation. To show the invertibility of L 2 , we claim that (Jm(Z/ 2 ))" L = {0} . 
Suppose if) G (/m(Z/ 2 ))" L . Then it is easy to check that L 2 [ip] = 0, -0 G X Q fl S . By Lemma 
12.11 we just need to show that i/j is bounded. To this end, we note that since ip G So, we have 
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ip(0) = e^^- 2 h ip(0) and k > 2 which may give ^(0) = . Let <j) = ip e^*. Then G X a and 
L 2 [(p] = 0. Hence G F Q and 

A0(z) = h{z)/ (1 + |^| 2+Q ) 1/2 ,V 2 6i 2 , 

where h G L 2 (R 2 ). Consequently, by Lemma 12.21 and Riesz representation formula, we may 
obtain 

\<t>{z)\ < Clog(l + |*|), Vzel 2 , 

and 

\d 6 (j){z)\ < CTog(l + |z|) , V^Gl 2 . 

Thus 

1^(^)1 , |£W>(*)| <CTog(l + |z|), Vzel 2 . (2.12) 
The equation I^V'] = 0, ^ = + \/— 1^2 can De written as 

d 2 , 2d. 



2 V^! - - 
7 

d 2 , 2d 



A^i + (1 - - - — «9^ 2 = 0, (2.13) 



A^ 2 + (1 - S 2 )^ 2 - ^ 2 + ^d e if>i = 0. (2.14) 
From (12121) and (12TT31 . Vi satisfies 

|A^i + (1 - 352)^1 - ^i| < ^ logr for r = |z| > 1 
Hence by comparison principle, 



d 2 . C 



^2 i^<2 



\Uz)\<j^-y (2-15) 
Here we have used the fact that Sa{r) ~ 1 as r — > oo. Similarly, we may have 

19^)1 < T ^-. (2.16) 

It remains to show that ^2 is bounded. Now we can use (I2.14p and (12 . 1 6[) to get 

d 2 C 

|A^ 2 + (1 - S 2 )^ 2 - -H < m - v for r = \z\ > 1 . 
In fact, ^2 satisfies 

d 2 2d 

-A^ 2 = /(*), where f(z) = (1 - S 2 - -)^ 2 + 
Since ip G E , we deduce that f(z) G X Q PI S and L 2 /(z) = 0. Since ^2(0) = 0, we have 

M*)=M0) + il- f \og r ^ 1 f(T)dr = ±- [ \ogP^-f(r)dr (2.17) 
2vrJ M 2 \z-t\ 2tt J r2 \z-r\ 

from which we obtain that 

<C\\f\\ Xa < +00 for zGR 2 . (2.18) 

Therefore by (I2.15P and (I2.18P , we may complete the proof of Lemma 12.31 □ 
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To study the vortex solutions, we perform the following key transformation: 

v = « d (z)e^W , Vd ( z ) = Sd (r) e^™ , 



(2.19) 



where ip = ipi + a/— 1^2 and ^ G 1R , j = 1,2. (Here we have assumed that the vortex occurs 
at the center.) Now we define 



Si[u, tp] := Au - u + u 3 + PuS 2 e~ 2 ^ 



S 2 [u,iP] := A^ + ^^- 



-is3(i - e- 5 ^; 



1|V# 



-l(3u 2 



We also write 
where 



S 2 [u,1>] = S 2 [u,0] + L 2 [iJ>] + N[u,4] 



Aip + ^^Vtfj - 2 v /r LS> 2 



(2.20) 



(2.21) 



A^i + 2 
+v /z T 



S d (\z\)\z\ 



Vipt - dve ■ v^ 2 



A^ 2 - 2£> 2 + 2 ( §^ ) V^ 2 + dVO ■ 



(2.22) 

V] = v^lV^I 2 - V=IS%(1 - e" 2 ^ - 2^ 2 ) , 
where V# = -(— sin cos 0). Then it is easy to see that solving (jl.5p is equivalent to solving 

S , i[m,V]=0, S , 2 [u,V] = 0. (2.23) 

Let ip = ipi + a/— 1^2, h = hi + y/—lh 2 . We may define two norms as follows: 



sup[|^i| + (1 + M)|v^| + (1 + M) 1+Q |^ 2 | + (1 + |z|) 1+t W 2 |] 



(2.24) 



and 



\h\\„ = sup[(l + |z|) 2+a (|/n| + \h 2 \ 

zeM 2 



Then we may show the following key lemma 
LEMMA 2.4. For any h G X a H So with < +00, i/iere exists a unique ?/> G F a fl £ 

L 2 [^] = ft (2.25) 

Furthermore, we have 

< C\\h\\„ (2.26) 
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Proof: Let h = — \/— 1 vjn and = — \J — 1 . Then ( 12.251) is equivalent to 

L 2 [$\=h (2.27) 

where /i satisfies h G Eo and 

ll^ll** < < +oo . (2.28) 

Consequently, (12.281) implies that L 2 |^| 2 (1 + |^| 2+a ) < +oo and hence G X a nT,Q. By Lemma 
12.31 (12.271) has a unique solution 0. Furthermore, as for (12.17j) . we obtain 



<f>{z) = — / log^^/(r)rfr (2.29) 

where |/(z)| (1+|2|) 2+Q < C||/i||**. From ([129]), we may deduce that |</>|+(l+|z|)| V0| < C\\h\\„ 
which implies that |^| + (1 + M)|V^| < C\\h\\„ i.e. |^-| + (1 + \z\)\Vip-\ < C\\h\\**,j = 1,2. 
To obtain better estimates for ip2, we may use the equation for ip2 and (12.221) . Then 

Mh - 2Sfah + O (jjj = O {\\h\U (1 + \Ar 2 ~ a ) , 

which gives 

<-A^ 2 + 2SU 2 < 



Hence we may use a barrier and elliptic estimates to get 

m\ + iv^l < 



Here we have used the fact that S^r) ~ 1 as r — > oo. Therefore we may complete the proof of 
Lemma EH □ 

3 Approximate Solutions and Error Estimates 

In this section, we introduce some approximate solutions and derive some useful estimates. 



Si[u,ij)] 
S 2 [u,ip] 



Let 

S[u,ip) = 
Let 

€j = le^^ '^, J = 1, k, Wj(z) := w(z - £,■), 

k 

Then we have 
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LEMMA 3.1. Fori large enough, we have 

||5i[^,0]|| L2(R2) + ||5 2 K,0]|| M <C(|/3|/ 2+a + e- 2Zsin f). (3.1) 

Proof, it is easy to check that 

Si[ui, 0] =Aui - ui + uf + f3ui\v d \ 2 



Due to 



we may obtain 



i^j 3 

/ < Ce" 21 ^' 1 < Ce- 4 ' sin f , Vi 7^ j , (3.2) 



||£ih,0]|| i2 < [e" 2 '^ + . (3.3) 

Here we have used the fact that £j's are vertices of regular /c-polygon with side length 21 sin ^ . 
On the other hand, we have 

k 

2|*-&| 



\S 2 [ Ul ,0]\ = \P\u*<cY,\P\e 

3=1 



and so by ffl~2Tj) 



||5 2 [ Ul ,0]||„<C|^|8up £ |^| 2+a e- 2 l^l <C|/3|^|0| 2+Q <C|/3||/| 2+Q . (3.4) 
Therefore by (13. 3p and (13. 4p . we may obtain (13.11) and complete the proof of Lemma [3.11 □ 

4 Liapunov-Schmidt reduction process 

Let X := (L 2 (M 2 ) x X a ) n E, F := (# 2 (M 2 ) x F a ) n S and L := ( - 1 I : Y -»• X , where 



L 



2 



la[0] = A0 — + 3m 2 and L 2 is defined by (12.211) . To solve (I2.23p . we first consider the 
following linear problem: Given / G L 2 (IR 2 ) R S 1; find (0, c) such that 

Z ^ = / + C ?T' 0e^ 2 (M 2 )ns l7 

where Si = {0 = 0(z) G R : (0, 0) G £}. 
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LEMMA 4.1. For each f G L 2 (R 2 ) n E 1; there exists a unique pair (0, c) satisfying H4-l\) 
such that 

\\4 H2 <C\\f\\ L2 . (4.2) 

PROOF. The existence and uniqueness may follow from Fredholm's alternatives. Now we 
prove (14. 2p by contradiction. Suppose H4.2[) is not true. Then there exist /3 n , l n G M, 4> n G 
# 2 (M 2 ) n Ei , /„ G L 2 (M 2 ) n Ei and c n such that (3 n -> 0, Z n -»■ +oo , 



L 2 



0, ||0, 



I// 2 



as n — > oo, 



and 



Ll(j) n — fn + C n , 



9/ 



0. 



(4.3) 

(4.4) 
(4.5) 



Let 4> n (z) := <f) n (z + £i). Then by (14.41) . <p n satisfies 
+3 



fe / k 

2 J] + 6 - 0) + ( ^ + & ~ &) 



J"=2 

/n(* + fl) +Cn 



</>n(z) 



(4.6) 



i=2 



We may multiply ( 14.4ft by -7^- and integrate over the whole space M 2 . Then by ( 14. 3p . it is 
obvious that 



dui 
dl 



(A(j) n - 4> n + 3u 2 4> n ) 



dui 
dl 



fn 



n— >+oo 



"ST ~ ~df l ~dl 
-> 0. 



So c n — > as n — > +oo. Hence by ( 14. 6p . as n — > oo, 



which satisfies A0 O — 0o+3u> 2 0o = 



in M 2 . Thus by [2n], 0o = % ^ — for some constants a\ and a 2 . Moreover, due to (j) n {z) = 

■ i OZj 

~ ~ ~ ~ <9u/ 

n (z), we have 0o(z) = 4>o( z )- Consequently, ai = and 0o = ai— — . On the other hand, by 

OZ\ 
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4>n{z) 



dzi 



n— »+oo 



3=2 

2 

<9zi 



Therefore ai = 0, 0o = 0, and <p n — > in L 2 oc (IR 2 ). Then we have 3uf(j) n — > in L 2 and hence 
||0n||# 2 — * which may give a contradiction and complete the proof. □ 

From Lemma 12.41 and 14. lj, we may obtain that 



Lemma 4.2. For 



h 

f2 



G X, there exists a unique 



in 



, c 6 7x1, such that 







fl 
h 



dui 

c | di 
o 



dui 



(4.7) 



Moreover, we have \\<p\\m + \\4>\\* < C(||/i||l2 + ||/ 2 ||**)- 

We may denote A as the inverse operator for Lemma \4. 21 i.e. A 
we have 

LEMMA 4.3. For I large satisfying $1.21\) . there exists a unique 



fi 
h 



in 







. Finally, 



G Y such that 



dui 

S[u l + (j) h ^ l ] = c{l)\ dl 





(4.8) 



Furthermore, 



*>i\\H2 + \m\\* 



Proof. This may follow from standard contraction mapping principle. We choose (0, ip) G 
B, where 



B 



eY:U\\ H ,+ 



<C(t 



-21 sin - 



\P\l 2+a ) 



and then expand 



S, [ Ul + M] = St [ Ul , 0] + U W\ + Ni [0, V] , 
+ 0, V'] = S 2 [u h 0] + L 2 [V] + N 2 [(f), 0], 
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where 

iV a [0, j>] = 3 Mi 2 + 3 + /? [0 ^ 2 e" 2 ^ + Mi S 2 (e" 2 ^ - 1)] , (4.9) 

and 

N 2 [</>,if>] = -V=Tp(2ui<f> + 2 ) + v 7 ^ |Wf - >/=TS2(l - e~ 2 ^ - 2tp 2 ) . (4.10) 
By (TOT)) and (0, ^) G £, we have 



^ + IHI*<|/3| 1 -' T , (0,— 



l 

100 



as (3 > sufficiently small i.e. I sufficiently large. Moreover, by (j4.9j) . (14.101) and the norms 
defined at (I2.24p . we see that 

||iVi|U 2( M 2 )<C|/3|||0|| L2 + ||0|| 2 2 + |/3|, (4.11) 

and 

\\N 2 \U*<CU\\ 2 + C\P\. (4.12) 

Now we can write (14.81) as 



S 1 [u l ,0]+N 1 [ ( p,ij} 
V> J " V 5 2 [«i,O]+JV 2 [0,V] 



.4 



Then as for the proof of Proposition 1 in [TB] , we may use a contraction mapping argument to 
obtain the desired result. Here we also need to use (14. lip . (14.121) . Lemma [3.11 14.11 and 14.21 □ 

5 Finding zero of c(Z) 

To prove Theorem it is enough to find a zero of c(l) in (14.81) . We multiply the first 

equation of (14. 8 p by — — and integrate it over the whole space M?. Then we obtain 

ol 

c{l) 2 = £ 2 [A(«, + 0,) - («, + 0z) + («, + 0O 3 ] ^ (5.1) 

+ /? / K + 0z)|^ + ^| 2 ^ 
= : h + I 2 , (5.2) 
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where 



h=P 
=P 

=P 
=P 



Note that 



2 dui 



Ui\v d \ A -^ + P / O(|0i| + K + 0i||Vj|) 
n ^^ + ( |/3|2/2+a + |/3|e ~ 2 ' Sin; 



01 



l-So+o(±)) ^ + 0(\(3\H^ + \(3\e-' 



2r 2 



(9/ 



dm 



u r-^ + o[\p\i- A + \p\ 2 i 2+a + w 



-2«sin£ 



dui 
dl 



E^E 

i =1 i^j 



dWj 
dl 



--0[e- 2lsin % 



P / u t 



d 2 dui 
r 2 dl 



-P / «i 



d 2 dui 
r 2 dl 



--P 



d 2 dw 
z + ^i| 2 dzi 



+ e 



-21 sin ; 



= -c 2 (3r 3 / w 2 + o(\p\r 5 + \p\ e 

where c 2 is a positive constant independent of j3 and I. So 
J 2 = -c 2 pl~ 3 I w 

For J 1; we have 



-2Zsin£ 



1 + 0(\(3\r 4 + \(3\ 2 l 2+a + \(3\e- 2lsin ^ 



h= I (A Ul -u l + uf)^+ f (A^-^ + 3n 2 0/)^ + O(|/5| 2 Z 2+a + |/5|e- 2isin f 



-E 



dui 
dl 



A 



<9u; dui 

~dT~~dT 



3u] 



dui 
dl 



^ + 0{\P\ 2 l 2+a + \P\e 



6k / w\w2 



^\ +0 {\P\ 2 l 2+a + \P\e- 21 ^ 

Z\ — I 



--6k w 2 (\ z -^\) w '(\ z -^\) w (\ z -^\)^^-dz + (\p\*P** + \p\e 

= -6fc / W 2 (| Z |X(H) W (^ + ei-6l)^r^ +0f|/5| 2 / 2+a +|/5|e- 2Zsin ^ 
-co •^(|6-6l) + 0f|/5| 2 / 2+Q + |/3| 



-2£ sin ; 



-21 sin? 



c (2sin^j .l-y 2 e- 2lsi ^ +0 (\P\ 2 l 2+a + \P\e- 2lsi ^y 
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where Cq is a positive constant. Here we have used the fact that (11.71) holds. 
Combining (15.51) and (15.61) . we see that 

c(i) « Cl r 1/2 e~ 21 sin i - c 2 (3r 3 + o (pr 5 ) . 

Moreover, we may choose lp such that 

f-VV^f =/3f- 3 . (5.7) 

Then 

^ = 2^nT| l0 4 + C3l0gl0 4 + C4 - 

Now we want to choose / such that I E (lp — 7, //? + 7) and c(l) = for some suitable 7 > 0. It 
is remarkable that 

c(lp - 7) « dft - ^-^e-^^f +^-5 - c^fe - 7 )" 3 + o(/? (log i)~ 5 ) > , 

if 7 is large enough. Similarly, c(jp + 7) < if 7 is large enough. Since c(l) is continuous 
in Z, then by the mean-value theorem, there exists lp E (lp — 7, lp + 7) such that c(/^) = 0. 

Consequently, the function ( Z/3 v /=^!f J =: ( J 9 J is a solution of (II. 5p . Furthermore, it 



w d e v ^ J \vp 

is easy to check that (up, Vp) satisfies all the properties of Theorem 11.11 Therefore we may 
complete the proof of Theorem 11.11 



6 Proof of Theorem 1.2 



We first consider problem ( I1.16P on a ball Br. 

Au — u + u 3 + (3uS 2 = 0, u = u(r),r < R 
AS-^S + S(l - S 2 ) + f3u 2 S = 0,S = S(r),r < R 
S(0) = 0,S(R) = l,u(R) =0,w>0,0<S<l 



(6.1] 



Our idea is to find a solution of (16. ip . and then let R — »■ +00. To this end, we consider the 
associated energy functional 



E R [u,S] 



\Vu\ 2 + u 2 ^ 



Br 



Br 



:i-s 



2\2 







Br 



WS\ 



2„.2 



d 2 



;S 2 



S z u 



Br 



U 



Br 



for u E H^(Br) and S E I R = {S E H\B R ) : S(z) = S{\z\), S(0) = 0, S(R) = 1}. 
Let the Nehari manifold be 



Af={(u,S)EH 1 (B R )xI R ,u>0,u^0: [ (\Vu\ 2 + u 2 )=! ( 
I Jb r Jb r 



(u 4 + (3u 2 S 2 ) 



(6.2) 



Then we consider the following energy minimization problem 



Cr = , ^ f ^rK 5 ] ' ( 6 - 3 ) 



and we have 



LEMMA 6.1. If (3 < 0, then cr is obtained by some radially symmetric function (ur,Sr). 
Furthermore, u R (r) < and S R (r) > for r > 0. 

Proof: We follow the proof of (1) Theorem 3.3 of [TH]. To this end, we define another energy 
functional 

E' R [u,S] = -j {\Vu\ 2 + u 2 ) + \ f (\VS\ 2 + ^S 2 ) (6.4) 
4 J b r 1 Jb r r 



R 



and another solution manifold 



W = \{u,S) e Hl(B R ) x I R ,u > 0,u^ : / (|Vw| 2 + u 2 ) < [ (w 4 + (3u 2 S 2 ) X . 
I Jb r Jb r J 

We consider another minimization problem: 

c' R = inf E'[u,S}. (6.5) 

Certainly, we have 

Cr < Cr. (6.6) 

Let (u n , S n ) be a minimizing sequence of c' R on A/"'. Replacing S n by min(S n , 1), we may 
assume that S n < 1. We may denote u* and S 1 * as the Schwartz symmetrization of u n and S n , 
respectively. Then (1 — S n )* = 1 — S*. By Theorem 3.4 of [TO] . 



(i-w< / (i-3D*«r (6-7) 

and hence due to j3 < 0, 

«) 2 (S;f<-p[ u 2 n S 2 n (6.8) 



On the other hand, we also have 



1 1 d 2 1 C 1 1 d 2 1 



W + ~ 4 (<) } = + 2^ ~ 4 <} 

[ (|v<| 2 + |vs;| 2 )< / (|vu n | 2 + |wg 2 ). 

Jb r Jbr 
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Hence we obtain 

E R K,S* n }<E R [u n ,S n ] 

and 

/ (|V<| 2 + «) 2 )</ («) 4 + /9«) 2 (^) 2 ). (6.9) 

J B R J B R 

Thus we may replace (u n ,S n ) by its symmetrization (u*,^*). Since S n < 1 and ^(Br) is a 
compact embedding to L A {B R ), we see that {u n S n ) — > (u^, Sr) weakly in if 1 (.Br) and strongly 
in L 4 (B R ), where (w^, S R ) G A/"' attains c R . So c R is attained. If (u R , S R ) G (A/"') — the interior 
of A/ 7 , then (w^, ,Sr) is a local minimizer of E' R and hence we have 

Au R — u R + (3u R S R = 0, V0 < r < R, and u R (R) = , 

which implies = since (3 < 0. This is impossible since from ( 16.91) and Sobolev embedding, 
we infer that f B u r > C > 0. Therefore (u R , S R ) G d(J\f') = Af and hence 

cr < E R [u R , S r ] = E R [u R , Sr] = c R . (6.10) 

Combining (I6.6P and (16.101) . we conclude that c R is attained by (u R , S R ). Then we have the 
following equality: 

Gr[ur,S r ]= [ (\Vu R \ 2 + u R -f3u R S 2 R -u R ) = 0. (6.11) 
Jb r 

Hence there exists a Lagrange multiplier X R such that 

VE r + XrVG r = 0. (6.12) 
Acting (I6.12p on (u R ,0), we may obtain 

I (\Vu R \ 2 + u 2 R - Pu R S R - u R ) + 2\ R j (\Vu R \ 2 + u 2 R - Pu 2 R S 2 R - 2u R ) = , 

J Bn J Bo 



and hence by (16.111) . 



-Ar / u% = , i.e. A 



R I U R - u > L-* 1 - A R 



□ 



Therefore, we may complete the proof of Lemma 16.11 

Theorem 11.21 is proved by the following lemma 

LEMMA 6.2. As R —> +00, (u R , S R ) — > (woo, S^) and (mqo, ^oo) is a solution of 111. 16]) . 

Proof: Since Sr < 1, Sr(0) = 0, we first show that ur is uniformly bounded, independent of 
R > 1. Actually, it is sufficient to show that J Br {\^ u r\ 2 + u r) < C, where C is a positive 
constant independent of R > 1. Let 

( l\VS\ 2 + f^S 2 + -(l-S 2 ) 2 



b r 2r 2 4 
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It is remarkable that GL Br may come from the conventional Ginzburg-Landau functional. 

From [2Tj . we may set Sr as the unique minimizer of GLr. Then for any u G Hq(Br), there 
exists t R such that (y/t^u, Sr) G TV, where t B is simply given by 



t R 



J BR (\Wu\ 2 + u"-f3Sl 



u 



(6.13) 



Thus by Lemma 16.11 and < Sr < 1 , 



cr < E R [ v / T R 'u,S R } = GL Br (S r ) + 



(X 



< GL Br (Sr) 

for all m G Hq(Br). Consequently, due to (5 < 0, 

Q? < GLb r {Sr) + Co 



LJv«l 2 + (i-/5K 



(6.14) 



where Co is a positive constant independent of R > 1. Here we have used the fact that 

r . - f BR \Vu\*+(l-PW 
hm mt — — 



(L R ^ 



< OO . 



By standard theory of Ginzburg-Landau equation (cf. [I]), we have 

GL Br [Sr] > GL Br [Sr\ . 
Combining (16.141) and (16.151) . we see that 

\ ! (|Vn R | 2 + 4-/55 2 4)-i f 

Z J Br, 4 JB, 



(6.15) 



<<C 3 



'Br 



and hence by the equation of ur, 



(\vu r \ 2 + 4) < c . 



(6.16) 



from which standard elliptic regularity theory gives that Ur < C. Thus we may obtain that 
(ur, Sr) — > (uqo, S^) which solves Au ro — !i M + + fiu^S^ = 0. Note that ur(0) > 1 and 
hence ^ 0. By the Maximum Principle, « M > 0. Similarly, < S^r) < 1 for r > 0. 
Therefore we may complete the proof of Theorem 11.21 

□ 
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7 Proof of Theorem 11.3 



In this section, we want to construct S 2 - valued map to get half-skyrmions by spike- vortex 
solutions obtained in Theorem 11.11 and 11.21 For simplicity, we firstly use spike-vortex solu- 
tions in Theorem 11.21 to find half-skyrmions. Let (u, v) be the radial spike- vortex solution in 
Theorem II .21 We may define a S 2 - valued map by 

1 ( Vl \ ( cos(0(r)) cos(d9) 

It = 2 ^ = cos(0(r)) sin(d0) (7.1) 

V / U 2 + V 2 + V i y u J y S in(0(r)) 

where v = V\ + y/^lvz, 



cos(0(r)) = ^L=, (7.2) 

sm(0(r)) = . (7.3) 

^u 2 + f 2 

Since both u and / are positive everywhere, the function is well-defined and single-valued. 
The map n can be decomposed into 

/ cos(dfl) \ / 

~n = cos(0(r)) sin(d#) | + sin(0(r)) ( 

V o 

Then it is easy to check that 

f d 

~n ■ (d x ~n A d y ~n) dx dy = / — -(f)'(r)cos((f)(r))dxdy 



and 



r 

l+oo 



-2?r(isin(0(r))|+~ = 2nd. 



i.e. the topological charge 



1 f d 

Q — — / ~n • (d x ~n A d y ~n) dxdy — - . (7.4) 
4vr J R 2 2 

Here we have used (17. 3D and the fact that w(0) > 0, u(+oo) = /(0) = and /(+oo) = 1. 

For the spike-vortex solution (u, v) in Theorem 11.11 since j3 > sufficiently small, the 
associated map n has the following form 

1 ( Vl \ ( cos (^) cos {dip) 

~™ = i o 9 =? 112 = cos ^) sin (^) I > ( 7 - 5 ) 
V u2 + V i+ V 2 \ u J \ sin(0) 

where <fi = 0(r, 9) and ip = ip(r, 9) satisfying 

cos(0M)) = ^4^' ( 7 - 6 ) 
\JU A + \v\ z 
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and tp = 9 + h, where h is a single-valued regular function satisfying h = 0{j3) as (3 — ► 0+. 
Here both u and \v\ may not have radial symmetry. Due to /3 > 0, we may apply the standard 
maximum principle on the first equation of the system (11.51) i.e. the equation of u. Then the 
solution u is positive everywhere so the function <fi is well-defined and single-valued. 

Now we want to calculate the topological charge Q as for (17.41) . By (17.51) . it is easy to check 
that 

n ■ (d x n A d y n ) = — r cos 0-1 — (06i/i r — 0r^e) cos • 
Hence by ( 17. 6p . (17. 7p and using integration by part, we may obtain 

1 f d 

Q = — ~n ■ (d x ~n A d y ~n) = - . 
4vr J R 2 2 

Here we have used the fact that u(0) > 0, u(oo) = 0, v (0) = and |u(oo)| = 1. Therefore we 
may complete the proof of Theorem 11.31 
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